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Abstract 

We construct a perfectly matched absorbing layer for stationary Schrodinger equa- 
tion with analytic slowly decaying potential in a periodic structure. We prove the 
unique solvability of the problem with perfectly matched layer of finite length and 
show that solution to this problem approximates a solution to the original problem 
with an error that exponentially tends to zero as the length of perfectly matched 
layer tends to infinity. 
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1 Introduction 



The Perfectly Matched Layer (PML) method, introduced in [1], is in common 
use for a numerical analysis of a wide class of problems. For some problems 
the convergence of the method has been proved mathematically, see e.g. [2-5]. 
In this paper, we introduce the PML method for the stationary Schrodinger 
equation in a "half-plane" with periodic boundary and Dirichlet boundary 
condition. We suppose that the potential q allows an analytic continuation to 
a cone on some distance from the boundary and q(z) uniformly tends to zero 
as z goes to infinity inside the cone. We include into consideration potentials 
decaying at infinity as slowly as z 11 , v < 0, or even 1/ Inz. Since the potentials 
are not compactly supported the modal analysis employed in [I] cannot be 
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used here, this leads to significant difficulties. Using the tools of complex 
scaling [8,9] we construct a PML of infinite length for the original problem 
supplied with some generalized radiation condition. The form of this radiation 
condition is similar to the pole condition PfTU] . The generalized radiation 
condition turns out to be equivalent to the classical radiation condition in 
the case of sufficiently rapid decay of the potential. As an approximation 
of a solution satisfying the original problem and the radiation condition, we 
take a solution to the problem with PML of finite length. We prove that the 
problem with PML of finite length is uniquely solvable and that the error of 
the approximation tends to zero with an exponential rate as the length of 
PML tends to infinity. The proof is based on weak statements of problems in 
weighted Sobolev spaces [7] and on a modification of the compound expansion 
method [6]. 

We consider the Dirichlet boundary condition as a boundary condition of the 
original problem and as an artificial boundary condition, however one can use 
the Neumann boundary condition instead. The approach is easily extended 
for this case, the results remain the same. 



2 Statement of the problem 



Let V be an upper "half-plane" in M 2 with smooth 27r-periodic boundary dV 
such that V C {(y,t) e R 2 : t > c} and dV C {(y,t) G R 2 : t < 0}. Let 
£ = {( y ,t) e V : \y\ < tt} be the periodicity cell of V. We set T ± = {(y,t) G 
V : y = ±7r} and T° = dS \ {T + U Y~}. As it usually is, the problem in V 
reduces to a quasi-periodic boundary value problem in the periodicity cell £, 
see Fig. 1. 



T 




Fig. 1. Geometry of the problem. 



We consider the stationary Schrodinger equation 



(A + k 2 + q(y, t))u(y, t) = F(y, t), (y, t) G £, 



(1) 



with the quasi-periodicity conditions 



aju(7T,*) 



e 2 ™diu(-ir,t), j = 0,1, (±7r,f)eT± 



(2) 
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and the Dirichlet boundary condition 



u(y,t) = 0, (y,t)er°. (3) 

Here a G [0, 1), d y = d/dy, and the parameter A; is a fixed real number that 
does not coincide with a threshold value, i.e. k 2 ^ (n + a) 2 for all n G Z. Let 
i£y denote the closed cone 

K$ = {zeC:z = T + e # t, < V < 0, * > 0}. 

We assume that the potential q in the equation (JTj) satisfies the conditions: (i) 
g is a bounded real-valued function in V, q(y, t) = q(y + 2n, t) for all (y, t) G V; 
(ii) for some T > and G (0, 7r/2) the potential q can be continued to an 
analytic in z (and 27r-periodic in y) function K x 9 (y, z) > q(y, z) G C, 
which uniformly tends to zero as \z\ — * +oo. 

We also make the following assumptions on the right hand side F of the 
equation ([1]): (i) F is in the space L l 2°(S) of locally square summable functions 
on £; (ii) for some T > and G (0, 7r/2) the function F is an analytic in 
^ function (— it, it) x 1^ E3 (y, z) \— > F(y, z) G C satisfying the uniform in 
ip G [0, 0] estimate 

/" + 00 /*7T 

/ / exp(2tr sin ij)\F (y, T + e i m 2 dydt < Const (4) 

JO J -TV 

with some r > 0. 



3 Radiation condition and complex scaling 



For all n G Z we set A^ = Tyk 2 — (n + a) 2 , where we take the main branch 
of the square root. Let 9t = {n G Z : |n + a| < |&|}. The finite set of points 
{A : A = A+ or A = A~, n G 91} consists of all points A^ lying on the real axis. 
The remaining points A„, n G Z \ 91 are on the imaginary axis. With every 
A^ we associate the function w^(y, t) = exp(i\^t + i(n + o)y). The functions 
satisfy the quasi-periodicity conditions (j2j) and the homogeneous equation 
([T]) with q = 0. If n G 91 then w+ is an incoming wave and w~ is an outgoing 
wave of the problem ([I])-© with q = 0. If n G Z \ 91 then is a growing 
mode and w~ is an evanescent mode of the unperturbed problem ([I])-©. Let 
be the angle for which the assumptions of Section [2] on the potential q and 
the right hand side F are satisfied. We introduce the open cone 

IC+ = {A G C : A = ip-e-^e (O,0),£ > 0} 

with the vertex i(3 G C and the angle 0, G (0, 7r/2). Denote by Hl(£), 
£ > 0, the weighted space with the norm ||eg-; where H l {£) is the 
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Sobolev space, and ep : (y, t) \— > exp /3t, (5 G R. We say that the parameter /3 
is admissible if (3 G (max{3A+ :nGZ\ 91}, 0) fl [— r sin 0, 0) and the cone /C^ 
contains all the points from the set {A+ : n G 91} , see Fig. 2. 



3A 




- incoming waves; 

- outgoing waves; 

- evanescent modes; 

- growing modes. 



Fig. 2. Arrangement of the cone )Ct for an admissible (5. 

Definition 1 A solution u to the problem ([I])-® satisfies the radiation con- 
dition if for some admissible j3 the solution u is in the space Hp(£) and the 
Fourier- Laplace transform u(y,\) = f^ 00 e~ tXt u(y,t) dt is an analytic in the 
cone K.0 function A i— > u(-, A) taking values in the Sobolev space H 2 (—tt, 7r). 

Theorem 2 Let the assumptions of Section^ be fulfilled, (i) If the homoge- 
neous problem (CO)-® has no nontrivial solution in the space H 2 (£) for any 
7 > then the problem ([I])-© has a unique solution satisfying the radiation 
condition, (ii) Let a solution u to the problem (EE])-® satisfy the radiation 
condition for an admissible (3. Then u satisfies the radiation condition for 
every admissible (3. 

We briefly discuss our formulation of radiation condition, for the details as well 
as for the proof of Theorem [2] we refer to [8] . One can note that our radiation 
condition looks quite similar to the pole condition introduced in [5|TU] as 
an equivalent and universal formulation of the classical radiation conditions 
for a wide class of problems. Formally, the only difference between the pole 
condition and our radiation condition is that we require the analyticity of the 
Fourier-Laplace transform in a cone instead of the half-plain. Nevertheless, 
the classical radiation conditions are not applicable under our assumptions 
on the decay of the potential. The introduced radiation condition should be 
considered as a generalization of the classical one. In the case of sufficiently 
rapid decay (say, with an exponential rate) of the potential q and of the 
right hand side F at infinity, our radiation condition is equivalent to the well 
known one: a solution satisfies the radiation condition if the principal term of 
its asymptotic at infinity is a linear combination of the outgoing waves. 

Let £ T = {(y,t) E £ : t < T}. By applying the complex scaling t — >• T + 
e l( ^(t — T) for t > T (complex change of variables) to the original problem 
(1-3), we obtain the problem 

(A + k 2 + q(y, t))v(y, t) = T{y, t), (y, t) G £ T , (5) 
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(aj + e^d 2 + qUy, t) + k 2 )v(y, t) = T{y, t), (y, t) E £ \ £ T , (6) 

^| x+ = e 2 ™^| x -, j = 0,1, (7) 

u = on T° (8) 

dfv(y,T-) = e-^div(y,T+), j = 0, 1, \y\ < vr, (9) 

where the equation (jUJ) is a jump condition, dfv(y,T—) and d(v(y,T+) denote 
the limits of dfv(y,t) as t tends to T from the left and from the right side 
correspondingly. The potential q% in (jSj) is defined by the equality 

q*(y, t) = q(y, T + e»(t - T)), (y, t)e£\W } (10) 

the right hand side T in ((SI),® is given by 

HvA = { FM ' M€£T _ (ID 

F(y,T + e i f>(t-T)),(y,t)eS\S T . 
The problem ©-© is elliptic for <p G (0, tt/2). 

Let i7^' a (£) denote the closed subspace in H}.{8) of all functions satisfying the 

quasi-periodicity condition w| T + = e 27rm M|x-- By H} f ,a (S) we denote the space 
of all functions u G H^ ,a (£) such that w = on T°. Then we introduce the 

space H~ l,a {£) as the dual space of H]l"{£) endowed with the natural norm 

\\T-H-^{S)\\ = sup{\(F,w) £ \ : w G H^(£), \\w;H^(£)\\ = 1}, (12) 

where (■, -)g is the extension of the inner product in L 2 (£) to the pairs (J 7 , w) G 
H~ l,a (E ) xifi'"(£). Let us note that if the right hand side F of the equation ([I]) 
satisfies the assumption (ii) from Section [2] for some i/;, T — T , and r = t 
then F satisfies the uniform in ip G [0, ip] estimate (HI) for every T > T , 
t < To, and ip < <p, see [9]. Without loss of generality we can assume that the 
potential q and the right hand side F satisfy the assumptions of Section [2] for 
some r > 0, G (0,7r/2), and for all sufficiently large positive T. Due to (JIJ) 
and (UK) we have T G H°(£) for all 7 < r sin 0, it is clear that H~ l ^{£) \\ < 
mH°{£)\\. 

Consider the variational statement of the problem find a function 

O 

v G H^' a (£) satisfying the equation 

—e~ l<t> J (d y v ■ dyW+d t v ■ d t w — (q{y, t) + k 2 )v ■ H)j dy dt 
- [ _(d y v ■ d y w+e~ 2i(f> d t v ■ d t w - (q£(y, t) + k 2 )v ■ w) dydt 

= e- i(t, {T, w) £T + (F, w) £ ^ w \/w G Hb%(£). 
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The variational form of the problem ©-Q generates the linear continuous 
operator 

H^iS) 3v^A lV efe H~ 1,a {£). (13) 

The proof of the following proposition can be found in [8]. 

Proposition 3 (i) Let the potential q satisfy the assumptions of Section^ 
and let T be a sufficiently large positive number. We define the potential 
for the equation ([6]) by the equality (jTUl) . If the homogeneous problem ([!])- 
([3]) has no nontrivial solution in the space H 2 {£) for any 7 > then the 
operator (|T3|) of the problem ([5]) -(Q yields an isomorphism if and only if 
| 7 | < min n6Z {3(e^A-)}. 

(ii) Assume that the potential q and the right hand side F of the problem ([T])- 
([3]) satisfy the assumptions of Section^ Let T be a sufficiently large positive 
number. We define the potential gT and the right hand side T of the prob- 
lem ©-([9]) by the equalities ffTUl) . ffTTj) . Let u be a (unique) solution to the 
problem ([I])-® satisfying the radiation conditions. Then a (unique) solution 
v G Hl ,a (E) to the problem is the analytic continuation of u in the 

sense that v = u on £ T and v(y, t) = u(y, T + e l ^{t — T)) for (y, t) G £ \ £ T . 



4 PML method. Rate of convergence and error estimate 

We search for an approximation in a domain £ L , < L < T, of a solution u to 
the problem ([I])-® subjected to the radiation condition. Since a solution v G 

Ho' a (£) to the problem ©-Q and u are coincident on £ T (see Proposition [31 
ii), one can search for an approximation of v instead of an approximation of 
u. The advantage is that v is in the space H^E), < 7 < min{S(e^A~)}, 
of functions "exponentially decaying" at infinity, while u ^ H l {£). It is clear 
that v has these properties because of the perfectly matched equation ([6]). In 
other words, the equation describes a PML of infinite length. 

We truncate the domain £ at a finite distance R > T. By T R we denote the 
boundary of truncation, = d£ R \ d£. Let us also set T ±,R = {(y, t) G T ± : 
t < R}. With the aim of approximating v by a solution v R to some problem 
in the bounded domain £ R , we introduce the problem 

(A + k 2 + q)v R (y, t) = F(y, t), (y, t) G £ T , (14) 

(fig + e -**flJ + fifa t ) + k 2 )v R (y, t) = T{y, t), (y, t) G £ R \ (15) 
d 1 y v R \ T+ ,R = e 2ma div R \ T -,R, j = 0,1, (16) 
V R = on T°, (17) 
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^ t v R (y,T-) = e-^div R (y,T+), j = 0, 1, (18) 
v R = G on T R , (19) 
where as an artificial boundary condition on we take the Dirichlet bound- 
ary condition. The equation ffT5l) describes a PML of the finite length R — T. 

Let H 1,a (S R ) denote the closed subspace in H 1 (S R ) of all functions satisfying 
the quasi-periodicity condition v\y+,R = e 27Tia v\ r ~,R and the boundary condi- 
tion v\y = 0. By H l,a (E R ) we denote the space of all functions v G H 1,a (S R ) 
such that v = on T° U T R . Then we introduce the space H~ 1,a (S R ) as the 
dual space of H l,a (£ R ). Consider the variational statement of the problem 
P11) - ([T51) : find a function v R G H 1,a (S R ) satisfying the equation 

-e _i * f (d y v R ■ d y w+d t v R ■ d t w - (q(y, t) + k 2 )v R ■ w) dy dt 

- f —(d y v R ■ d y w+e- 2i,t) d t v R ■ d t w - (q+(y, t) + k 2 )v R ■ w) dydt 

= e~ l \T, w) £T + {F, w) £RW G H 1 ^^) 

and the boundary condition v R — Q on T R . The variational statement gener- 
ates the linear continuous operator 

H l ' a {E R ) 3v R ^ A R v R = {JF, g} g H- l > a (£ R ) x H l/2 > a (T R ), (20) 

where H l l 2 ' a (T R ) is the space of traces on of the functions from H 1,a (£ R ). 

Proposition 4 Let T be a sufficiently large positive number and <fi G (0, vr/2). 
Assume that for all 7 > there is no nontrivial solution to the original ho- 
mogeneous problem (CO)-® in the space H 2 {£). Then there exists Rq > T 
such that for all R > Ro the problem (fl4l) -( [191 with right hand side {J 7 , Q} G 
H~ 1,a (E R ) x H 1 / 2 ' a (T R ) admits a unique variational solution v R G H 1,a (£ R ). 
The estimate 

||^ ; || < C{\\T ] H- 1 ' a {8 R )\\ + \\Q] H 1/2 ' a (T R )\\) (21) 

is valid, where the constant C does not depend on R > Ro. 

PROOF. The proof is carried out by a modification of the compound ex- 
pansion method |6j. In other words, we find an approximate solution to the 
problem (fl4^ - (fl9|) compounded of solutions to first and second limit problems. 
As the first limit problem we take the scaled problem (J5j) — (|9j) . The second limit 
problem is the elliptic problem with constant coefficients 

(d 2 + e- 2 ^d 2 + k 2 )l) 2 (y,t) = F 2 (y,t), (y,t) G IT* 

^U 2 (vr,t) = e 2 ™^U 2 (-7r,t), j = 0, 1, t < R, (22) 

U 2 = G 2 onT fl , 
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where II* = {(y,t) : y G (-n,n),t < R}. 

Let us define the functional spaces for the problem (122|) . By if* (II*) we denote 
the weighted Sobolev space with the norm if^I!*)!! = e~ 7 ||e 7 w; if^II*)!!. 
The space if*' a (n*) is the closed subspace in H^(H R ) of all elements satisfying 

u(n,t) = e 2ma u(—ir,t) for t < R. The space ifl' 7 (n*) consists of all elements 

u G H^ a (U R ) having the traces u\ t r = 0. We set if- 1 -^!!*) = (ili' 7 ( n -))*> 
the space ii^ 1,a (IT*) is provided with the natural norm, cf. (TT2"]) . Consider 
the variational statement of the problem ([22]) : /md a function U 2 G if 7 ' Q (II*) 
which satisfies the equation 

I (-d y U 2 -d y w-e- 2i ^d t U 2 -d t w+k 2 U 2 -w) dy dt = (F 2 ,w) n « Vw G ifi' 7 (II*) 

and £/ie boundary condition U 2 = G 2 on Y*. To the variational form of the 
problem (|22|) there corresponds the linear continuous operator 

H^(Tl R ) 9U 2 k A„ 7 U 2 = {F 2 , G 2 } G ifI 7 ' a (II*) x If ^(T*). (23) 

As is well known [TJ, the operator A_ 7 is Fredholm (i.e. the range of the 
operator is closed, kernel and cokernel are finite-dimensional) if and only if 
there are no numbers e^A^, n G Z, on the line {A G C : 3A = 7}. Suppose 
that 7 G [0,min n {3(e^A~)}). Then e^A± i {A G C : 3A = 7}, n G Z. 
The solutions to the homogeneous problem (|22|) are easily found in an explicit 
form, one can see that they do not belong to the space ifl 7 (n*). Analogously, 
we consider the formally adjoint to (|22|) homogeneous problem, and check that 
it has no solution in the space if 1 , (II*). Therefore, if 7 G [0, min^Q^e^A")}) 
then the operator (1231) implements an isomorphism, a variational solution U 2 G 
if_' 7 (n*) of the second limit problem (I2"2"j) with right hand side {F 2 ,G 2 } G 
ffI 7 ' Q (II*) x H x l^ a {T R ) satisfies the estimate 

||U 2 ;tf^( n -)H < C(||F 2 ;fC^(II*)|| + \\Q 2 -H l ^ a {T R )\\). (24) 

The constant C in (1241) is independent of R because the problem (122]) reduces 
to the same problem with R = by the shift t 1— > t + R, and the norms in (124]) 
are invariant with respect to ff, e.g. ||U 2 ; ffj(n*)|| = ||U 2 (- + R); ff 7 (II° )||. 

Now we are in position to construct the approximate solution. Let \ be a 
smooth cut-off function on the real line, x{t) = for £ > 1 and x{t) = 1 
for t < —1. We denote XR/2(t) = x(t ~ -R/2), t G R. For a sufficiently large 
i? > we set Fx = XR/2F and F 2 = (1 — XR/2)^ "> where (?/,£) G £* and 
JF g H~ 1,a (£ R ) is the right hand side of the problem (|T4^ - (IT51) . We extend 
the functional F x (the functional F 2 ) by zero to all t > R (to all t < 0). It is 
clear that T = Fi + F 2 , for all 7 G [0, min n {3(e^A~)}) we have 

HFijif- 1 ' ^)!! < e^'^W^H^^W, (25) 
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||F 2 ;#l}' a (II*)|| <e^ R/2+1) \\r ] H- 1 ^{8 R )\\. (26) 

Let Ui G H} / ' a (S) be a (unique) solution to the first limit problem (j5j)-((9]) 
with the right hand side Fx, and let U2 G if^' a (Il^) be a (unique) solution 
to the second limit problem f[2"2"j) with the right hand side F 2 and G2 = 0, 
where Q is the same as in (18), g G H 1/2 ' a (T R ). Due to the first assertion of 
Proposition [31 the estimate 

||Ui;^(5)|| ^C^H^^W (27) 

is valid. We define the approximate variational solution Y G H l,a (S R ) to the 
problem ffT4^ - ffl~9]) by the equality 

Y(y, t) = X (t - 2R/3)[i 1 (y, t) + (1 - X (t ~ i?/3))U 2 (y, t), (y, t) G £ R . 

By setting 7 = in the estimates ([24]), §ZS§, ([26]), and (J2ZD, we derive 

Hr;^ 1 '"^)!! < Constdl^if- 1 '"^)!! + ||£;# 1/2 ' Q (r R )||) (28) 

with some constant independent of J 7 , Q, and R. 

On the next step we estimate the discrepancy that Y leaves in the right hand 
side of the problem (|T4 l -( [T9|) . in other words, we estimate the value 

\\A R Y -{r,g} ] H- 1 ^{S R ) x £T 1 / 2 ' a (T fl )||; 

here A R is the operator ( J20l . Recall that by Ay and v4_ 7 we denote the oper- 
ators ( fT3l) and (|23l) of the first and second limit problems. It is clear that the 
mappings 

Hl' a {E) 3 Ui 1- AX2/3/JU! G H ^ ,a (£ R ), 
H l > a (II R ) 3 U 2 ^ A (l - Xi?/ 3 )U 2 G iT 1 '^*) x H l '^{T R ) 

are continuous. We have 

AY = {^oX2fl/ 3 Ui,0}+{g^(l-Xii/3)U2,0}+^o(l- XB/ 3)U 2 = (29) 

where we assume that the function g^(l — XR/3) is extended to £ R by zero. 
From §M) it follows that 

{^, Q} - {T, Q} = {Fx + [A, Xafl/s]Ui, 0} + {<#(1 - Xr/3 )\J 2 , 0} 

+ {F 2 ,g}-[A , X R/3]^2-{J r ,G} (30) 
= {[A,X2fl/ 3 ]Ui,0} + {4(1 - x«/ 3 )U 2 ,0} - [A ,Xi?/ 3 ]U 2 ; 

here [a, 6] = ab — ba. The term [y4 , X2i?,/ 3 ]Ui is equal to zero outside of the set 
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{{y,t):ye [-7T, n],te [2R/3 - 1, 2#/3 + 1]}. We get 

\\[A , X 2R/3W,H-^(£ R )\\ < Ce- 2 ^/ 3 ||U i; ^(£)|| 

< Ce-^ R/& \\T-H- x ^{S R )\l 

where the last estimate is a consequence of the estimates fl25|) and ff27j) . the 
constant C does not depend on R. A similar reasoning together with (T24l and 
( 1261) leads to the estimates 

||([^o,Xi ? /3]U 2 ) 1 ;^- 1 ' Q (^)|| < Ce- 2 ^/ 3 ||U 2 ;^( n -)H 

< Ce- 2 ^/ 3 (||F 2 ;^' a (n fi )|| + tf 1 / 2 ^) ||) (32) 
< Ce-^dl^; H~ l > a (£ R )\\ + ||0; H l ' 2 > a (? R ) ||) 

for the first component of the pair [A , = {([A), Xr/3]U 2 )i, 0}. At last, 

due to our assumptions on the potential q (see Section [2]), we have 

||4(1 - xr/s) U 2 ; H~^(£ R ) || < c(R) || U 2 ; H^ a (U R ) || , (33) 

where c(i?) tends to zero as R — > +oo. From (|33|) and the estimates (1241) . (I26I) 
with 7 = 0, we see that 

\\4(i -xr/^h-^{£ r )\\ < Cc(r)(\\^,h-^(£ r )\\ + ||f?; jy^cr*)!!). 

(34) 

Taking into account the equalities (1301 and the estimates fl3Tl) . fl32|) . and fl34|) . 
we arrive at the estimate 

WF-F-H-^i^W < C{R)(\\F;H- x > a {8 R )\\ + \\G\ ff 1/2 ' a (T fl )||), (35) 

where C(-R) does not depend on {J 7 , 0} and C(R) — > as R — > +oo. 

We first assume that = 0. In this case we have JF = JF + D(i?).F with some 
operator 0(R) in H~ 1,a (£ R ), whose norm tends to zero as i? — > +oo. Hence for 
a sufficiently large Ro and for all R > Rq we have |||D(i2)|||ij < |||D(-Ro)|||i? < 1' 
where ||| • |||# stands for the operator norm in H~ 1)Ci (£ R ). There exists the 
operator (/ + D(R))^ 1 : H~ 1,a (£ R ) — > H~ 1,a (£ R ), the norm of this operator 
is bounded by the constant 1/(1 — |||D(i?o)|||_R ) uniformly in R, R > R . 
We set T' = (I + D(i?)) -1 JF. In the same way as before we construct the 
approximate solution Y to the problem ffT4l) - (fT9]) . where T is replaced by 
T' and = 0. Then A R Y = {.F,0}, the estimate (TJHD holds with = 0. 
This proves that in the case = the problem (fT^ - (fT9]) has a solution 
v R G H 1,a (£ R ) satisfying the estimate (l2~Tj) . In the case Q ^ we find an 
exact solution v R to the problem (TT^l) fTITJT) in the form Yi — Y 2 . Here Y"i is 
the approximation solution of the problem (fl4|) - (fl9|) with the right hand side 
{.F, 0}, and F 2 is the approximation solution of the problem (TT4l) (n~9|) with 
= and T replaced by (I + Q{R))-\f -J 7 ). We have A R Y X = {? ', 0} and 
A R Y 2 = {T-T, 0}. Now we see that in the case ^ the problem (THD-CT 
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also has a solution v R G H 1,a (S R ) satisfying the estimate (12TT) . Indeed, by the 
proved case the estimate (12TT) is valid for v R = Y 2 , Q = 0, and T replaced 
by T — T . This together with the estimate (I3"oT) . and the estimate (1251) for 
Y = Yi, leads to (12T|) . To prove the uniqueness of the solution t> fi it suffices 
to apply the same argument to the formally adjoint problem. Proposition H] is 
proved. 

Theorem 5 Let T be a sufficiently large positive number. Assume that the po- 
tential q and the right hand side F satisfy all the assumptions of Section^ the 
homogeneous problem ([I])-© has no nontrivial solution in the space H^(S), 
e > 0. We define the potential q^ in ( fl5l) and the right hand side T of the 
equations (1T41) . (fl5l) by the equalities (TIP]) and (TTTj) . Let u denote a solution 
to the original problem (JXJ) — (J3|) with radiation conditions. Then there exists 
Rq > T such that for R > Rq a (unique) variational solution v R to the prob- 
lem (fT4l) - ffl9|) . where Q = ; converges to the solution u in the domain S T in 
the following sense 

Wu-v^H^S^W^Ce-^W^H^l R>R , (36) 

where the constant C is independent of R and T , and 

7 G (O,rsin0] n (0,min{3(e^A;)}) (37) 

with the same r > as in (TJJ . 

PROOF. Due to Proposition [3] a unique solution u to the problem (CD)-® 

with radiation conditions and a unique variational solution v G Hl ,a {8) to the 
scaled problem ©-([9]) are coincident on E T . Thus the estimate fl36l) is valid if 
and only if it is valid with u replaced by v. The difference v — v R G H 1,a (S R ) 
satisfies the problem fTHl) - f[Tgi) with the right hand side T = and Q = v\ t r. 
It is clear that 

IMxh;^ 1 / 2 ' ^)!! < e^ R \\v;H^(£)\\. 

This together with the first assertion of Proposition [3] and the assumption (Tj0) 
leads to the estimate 

Wg-H^IT^W < Ce-i R \\F;H°(E)\\ 

with the same restrictions on 7 as in ( 137|) . Then by Proposition H] we have 

\\v-v R ;H\S R )\\ < C\\g-H l ^ a {T R )\\ < Ce^ R \\T] H°{£)\\, R > R . 

Theorem [5] is proved. 
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Remark 6 It is quite possible that Theorem^ (as well as Propositions^ and 
0) remains valid without the assumption on the largeness of the parameter T. 
But we suppose it all the same because our proof of Proposition^ is essentially 
based on this assumption; see J2}/. 
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